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Abstract 

We perform the first study for the heavy-heavy-light quark (QQq) potential in SU(3) 
quenched lattice QCD with the Coulomb gauge. The calculations are done with the 
Ph, standard gauge and 0(a)-improved Wilson fermion action on the 16 4 lattice at 

j3 = 6.0. We calculate the energy of QQq systems as the function of the distance R 
between the two heavy quarks, and find that the QQq potential is well described 
£S) . with a Coulomb plus linear potential form up to the intermediate distance R < 0.8 

fm. Compared to the static three-quark case, the effective string tension between 
the heavy quarks is significantly reduced by the finite-mass valence quark effect. 
This reduction is considered to be a general property for baryons. 
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1 Introduction 



The inter-quark interaction is one of the fundamental and essential properties 
linking elementary physics and hadron physics. In particular, the three-quark 
interaction in baryons is characteristic and complicated, reflecting the non- 
trivial gluonic dynamics based on the SU(3) gauge symmetry. In addition, the 
three-quark system has also large varieties of the quark motion, configuration, 
and so on. The heavy-heavy-light quark (QQq) system is a suitable material 
to investigate such a light-quark effect on the three-quark interaction. 
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Fig. 1. (a)The gauge-invariant QQq Wilson loop. The wavy line represents the 
light-quark propagator and the straight line the heavy-quark trajectory. (b)The 
"wall-to-wall QQq Wilson loop." The gray wavy lines represent the wall-to-wall 
quark propagator, which propagates from the whole space at one time to that at 
another time. 

In 2002, the first doubly charmed baryon, H+(cfcc), was experimentally ob- 
served at SELEX, Fermilab PQ. In this experiment, a decay process — > 
A+K~it + was observed, and its mass was measured about 3519 MeV. In an- 
other experiment, a decay process — > pD + K~ was also confirmed [2]. 
Doubly charmed baryons are also theoretically investigated in lattice QCD [3] 
and other approaches [4ll5] . Since c quark is much heavier than d quark, this 
system can be idealized as the three-body system constructed with two static 
quarks and one finite-mass quark moving around. 

Motivated by these considerations, we investigate the QQq system in quenched 
lattice QCD. We extract the QQq potential VQQ q (R), which is defined as the 
energy of the QQq system in terms of the inter-heavy-quark distance R. In lat- 
tice QCD, the quark- ant iquark (QQ) potential [6] and the static three-quark 
(3Q) potential [7J are already found to be described as the linear confine- 
ment potential plus the one-gluon-exchange Coulomb potential. The confine- 
ment potential is proportional to the length of the gluonic flux connecting 
the quarks, and the string tension is about 0.89 GeV/fm. In contrast to these 
static quark potentials, our QQq potential includes not only the gluonic ef- 
fect but also the nontrivial finite- mass valence quark effect. Thus the QQq 
potential behavior would have some difference from the static cases. 



2 Formalism 



The QQq Wilson loop is denned in almost the same manner as the static 3Q 
Wilson loop [7J. The difference is that the light quark is treated as the quark 
propagator K~ x . The gauge- invariant QQq Wilson loop with the spatial size 
R and the temporal size T is defined as 



2 



Table 1 

The lattice parameter j3 = 2N c /g 2 , the corresponding lattice spacing a, the sweep 
numbers (A'therm; -^sep) of the thermalization and separation for updating the gauge 
fields, the smearing parameters (a,N smT ), and the clover coefficient c. 



a [fm] 


lattice size 


-^therm 


N scp 


a N smT c 


6.0 0.10 


16 4 


10000 


500 


2.3 40 1.479 



Table 2 

The correspondence between k and the used gauge configuration number N con f . The 
list shows the pion mass m^, the p meson mass m p , and the approximate constituent 
quark mass M q ~ m p /2. The meson masses are obtained from the meson correlator 
with the wall source and the point sink. The statistical error is estimated with the 
jackknife method. 



K 


-^conf 




m p a 


M q 


0.1200 


1000 


1.446(1) 


1.472(2) 


1.5 GeV 


0.1300 


300 


0.900(2) 


0.949(1) 


1 GeV 


0.1340 


300 


0.643(1) 


0.716(1) 


700 MeV 


0.1380 


1000 


0.304(1) 


0.467(2) 


500 MeV 



W QQq (R, T) = -e abc e def U l ad UllK- f \ (1) 

with the path-ordered product of the link variables U k = P exp(ig J r dx^A p ) 
along the heavy-quark trajectory Y k {k = I, II). The subscripts a,b,...,f are 
color indices. The schematic figure is shown in Fig. Q](a). The QQq potential 
is obtained as 



V QQq (R) = - lim i HW QQq (R, D>. (2) 

1 — ^oo ± 

The symbol ( ) means the expectation value integrated over the gauge field. 
All we have to do is to get the expectation value from lattice QCD for several 
values of R and to give a suitable function form of VQQ q (R). 



3 Simulation details 



We generate the SU(3) gauge configurations with f3 = 6.0 and 16 4 isotropic lat- 
tice at the quenched level. We adopt the standard plaquette gauge action, and 
the pseudo-heat-bath algorithm to update the gauge field. Periodic boundary 
conditions are imposed on the space-time boundaries. We apply the smearing 
method in Refs. [7]I8] to the spatial link variables of the QQq Wilson loop. The 
smearing method changes a stringy link to a spatially-extended flux tube, and 
enhances the ground-state component without changing the physical content. 
The method has two parameters, a real parameter a and the iteration num- 
ber N smr , and our choice of a and iV smi . is based on the static 3Q case [7]. 
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Fig. 2. The typical example of the effective mass plot: v(R = 8,T) with 
(Ri,Ru) = (4,4). The upper data correspond to the heaviest case, k = 0.1200, 
and the lower data the lightest = 0.1380. All the scales are measured in 

lattice unit. 

These simulation parameters are summarized in Table [H The lattice spacing 
a ~ 0.10 fm is determined so as to reproduce the string tension of the QQ 
potential to be 0.89 GeV/fm . We use the lattice unit for most part of the 
paper. 

For the light-quark propagator, we adopt the clover fermion action, which is 
the 0(a)-improved Wilson fermion action [9]. The clover coefficient c in this 
action is given from the mean field value uq of the link variable for the tad- 
pole improvement. We determine c and Uq from the ensemble average of the 
all plaquette values P^in) as c = 1/wq and u Q = (J2 n J2fj,>u |ReTrP Atl/ (n)) 1//4 . 
The measured value of u is 0.87779(2) in our case. To investigate the light- 
quark-mass dependence, we take different four light-quark hopping parame- 
ters, k = 0.1200, 0.1300, 0.1340, and 0.1380. Their correspondences to the 
light-quark masses are shown in Table [2j The constituent quark mass M q is 
roughly estimated with the half of the p meson mass. 

In Eq. ([T]) and Fig. UK a), the QQq Wilson loop is defined as a single gauge 
invariant loop. To reduce the statistical error, we adopt the following pre- 
scription. The light-quark propagator is spatially extended as the wall source 
and the wall sink. (The "wall" means the average over all spatial sites.) This 
propagator is averaged over all the spatial sites n src of the source and all the 
spatial sites n sin k of the sink at the fixed time separation T, as i^~ all (T) oc 
EnsrcEn 8fak ^ -1 ( n src,^sink,r). The "wall-to-wall QQq Wilson loop" is con- 
structed from this wall-to-wall quark propagator and the heavy-quark tra- 
jectories, and its schematic figure is depicted in Fig. [TJ(b). Because such a 
propagator is independent of the spatial position, we can easily sum up the 
parallel translated wall-to-wall QQq Wilson loops in the whole space. This 
summing up drastically suppresses the statistical error, owing to the large 
statistics, e.g., 16 3 times larger in our case. In the gauge invariant formalism, 
only gauge invariant components in wall-to-wall QQq Wilson loops remain. 
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Table 3 

The lattice QCD results for the Coulomb-gauge-fixed QQq potential Vqq q at ft = 
0.1380. R and (Ri,Ru) denote the loop size defined in Fig. [TJ The results with 
different fit ranges of T are also shown. All the values are in lattice unit, and the 
statistical error is estimated with the jackknife method. 



R 


Rji) 


V QQq inT=[4,8] 


V QQq inT=[5,8] 


1 


(0,1) 


0.877(2) 


0.873(2) 


2 


(0,2) 


0.971(7) 


0.959(9) 




(1,1) 


0.969(8) 


0.958(10) 


3 


(0,3) 


1.047(4) 


1.045(7) 




(1,2) 


1.045(4) 


1.043(8) 


4 


(0,4) 


1.083(11) 


1.067(17) 




(1,3) 


1.079(10) 


1.063(16) 




(2,2) 


1.078(10) 


1.063(15) 


5 


(0,5) 


1.136(6) 


1.122(3) 




(1,4) 


1.131(6) 


1.117(4) 




(2,3) 


1.130(6) 


1.116(5) 


6 


(0,6) 


1.170(13) 


1.151(24) 




(2,4) 


1.157(16) 


1.136(30) 




(3,3) 


1.157(16) 


1.136(31) 


7 


(0,7) 


1.219(21) 


1.220(50) 




(3,4) 


1.207(24) 


1.209(60) 


8 


(0,8) 


1.262(11) 


1.283(21) 




(4,4) 


1.255(6) 


1.271(10) 



When we calculate the QQq potential in the gauge invariant way, we find that 
the statistical and systematic errors are severely large, especially for lighter 
quark mass case. For the error reduction, we fix the gauge configurations with 
the Coulomb gauge. We should note that, in the wall-to-wall QQq Wilson 
loops with the Coulomb gauge, gauge variant components also remain due to 
the nonlocal nature of the gauge. However, the Coulomb gauge fixing empir- 
ically does not affect the long-range physics, such as the string tension |10j . 
In this paper, we are mainly interested in the long-range behavior of the QQq 
potential, and investigate it with the Coulomb gauge. We have also calculated 
in the Landau gauge, and obtained the same result as the Coulomb gauge. 



4 Results 



To decide the fit range of T, we define the effective mass 



v(R,T) = In 



(W QQq (R,T)) 
(W QQq (R,T + l)) 



(3) 
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and seek its plateau region against T. If the state is dominated by a single 
component, v(R,T) is independent of T. Typical cases are plotted in Fig. [2j 
All the statistical errors are estimated with the jackknife method here and 
below. We see that, in T > 3, the effective mass is approximately flat and 
thus the ground state component dominates. By fitting (Wqq q ) with a single 
exponential form Ce~ VQQqT , we obtain the QQq potential values Vqq q of dif- 
ferent loop sizes, as partially listed in Table [31 As mentioned above, Vqq q is 
almost unchanged for the different fit range of T. 

Remarkably, Vqq q depends only on R(= Ri + Ru) and not on the combination 
of (Ri,Rn). This property is suitable for the potential calculation for the 
following reason. If the fit range of T is large enough, the expectation value of 
the Wilson loop is dominated by the ground-state component, and does not 
depend on its condition of the source and sink. Then the QQq potential does 
not depend on (Rj, Rn) or more arbitrary choice of the junction points. 

To obtain the functional form of VQQ q (R), we consider the fit function, 



V QQq (R)=a eS R-^ + C eS , (4) 

as the analogy of the QQ potential. The subscript "eff" means the effective 
values including the light-quark effect. We find that this function is fairly 
suitable for Vqq q , at least in this calculated region of R, i.e., R < 0.8 fm. The 
best-fit parameters are shown in Table HJ and the corresponding potential form 
is depicted in Fig. [3j In the static QQ and 3Q potentials, the string tension 
and the Coulomb coefficient are obtained as 



a QQ - - 0-045, -Aqq ~ A 3 q ~ 0.13 (5) 

in the lattice unit at (3 = 6.0 [7j. The effective Coulomb coefficient A c q is 
almost the same value as A 3 q. As for the effective string tension a e g, in the 
heaviest valence quark case, a e s does not deviate from <t 3 q within the statis- 
tical error. However, a e g obtained in the lightest quark case is significantly 
reduced from <j 3 q, 



Ocff < cr 3Q , (6) 

provided that the finite volume effect is negligible. Therefore, when the valence 
quark is light, its effect reduces the effective string tension a eS from the string 
tension a 3 Q in the static 3Q case. 
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Fig. 3. The lattice QCD data of the Coulomb-gauge-fixed QQq potential Vqq q 
plotted against the inter- heavy-quark distance R for the different four values of k. 
The solid curves are the best-fit functions of Eq. All the scales are measured in 
lattice unit. 

5 Discussion and summary 



In our QQq potential, the significant result is the reduction of the effective 
string tension. To understand this result, we compare the string tension with 
the effective string tension. In the QQ case, the color flux-tube length corre- 
sponds to the distance between the quark and the antiquark, and thus ctqq is 
the proportionality coefficient of the inter-quark distance in the confinement 
potential. Similarly, in the 3Q case, a^Q is the proportionality coefficient of the 
flux-tube length that minimally connects the three quarks [7|lllll2j . However, 
in the QQq system, the inter-quark distance R is not the flux-tube length. The 
QQq flux-tube length is determined by complicated light-quark dynamics, and 
therefore VQQ q (R) does not have to be proportional to R. The actual lattice 
QCD results suggest that the long-range QQq potential is proportional to R at 
least for R < 0.8 fm, but the effective string tension o c s as its proportionality 
coefficient is reduced by the existence of the light quark, which is conjectured 
to originate from the difference between R and the flux-tube length. 

Now we consider the asymptotic behavior of the QQq potential. In the large R 
limit, the confinement potential between the two heavy quarks will dominate, 

Table 4 

The best-fit values of <7 e ff, Arff> and C c s in Eq. and their \ 2 over the degree of 
freedom A^of- All the values are in unit of the lattice spacing a ~ 0.10 fm. 



K 


C"cfT 




C c s 


X 2 /N doi 


0.1200 


0.045(2) 


0.12(2) 


1.49(2) 


1.31 


0.1300 


0.038(4) 


0.13(2) 


1.23(3) 


1.18 


0.1340 


0.037(4) 


0.13(2) 


1.12(2) 


1.11 


0.1380 


0.037(2) 


0.13(1) 


0.97(1) 


1.16 
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and the color flux-tube length approximately equals to R. Therefore a e g would 
approach to a^g, which will be confirmed with a larger- volume calculation. 
There is another asymptotic behavior we can check, that is, the limit in which 
the light quark is "heavy". If the light-quark mass goes to infinite, the QQq 
system should correspond to the static 3Q system. It means that a e s becomes 
larger in the larger light-quark mass and cr e g = <j^q in the infinite mass limit. 
We can confirm this behavior from the lattice result in the heaviest case, 
K = 0.1200, in Table H 

In our another work [T3], we investigate the same QQq potential in a non- 
relativistic quark model. It reproduces the lattice QCD result, and enables us 
to understand the reduction mechanism of a e s- We confirm that a geometrical 
difference between the flux-tube length and R is essential for the reduction of 
cr e flf, as conjectured above. 

We mention another possible light-quark effect on the interquark potential, 
i.e., the sea quark effect. The sea quark effect causes the string breaking, 
which is the disappearance of the string tension at a long distance. The sea 
quark effect is obtained with the unquenched lattice QCD calculation, and is 
found to be important at the long range where the flux-tube length is larger 
than about 1 fm [H] . The effective string tension in the realistic QQq system 
would be also affected by such a sea quark effect in further large- R region. 

Although our lattice calculation is restricted to the specific gauge, we come 
to the following general conclusion about the effective string tension, which 
would not be affected by the nonlocality of the Coulomb gauge. Owing to 
the existence of a mobile quark, the effective string tension a c g between the 
other two quarks can be reduced from the string tension in mesons <Jqq or 
baryons a^Q. Since the reason for the reduction is fairly simple and general, 
this argument holds for not only QQq systems but also three finite-mass quark 
systems, i.e., ordinary baryons, such as a nucleon. In addition, this can be 
also applied to the multi-quark system including light quarks [15] . For more 
quantitative calculation for realistic hadrons, we need careful consideration 
about realistic quark masses, finite- volume effects, sea quark effects, and more 
complicated valence quark effects. 

In summary, we have studied the QQq potential in SU(3) quenched lattice 
QCD. For the error reduction, we have fixed the gauge field with the Coulomb 
gauge, and investigated the long-distance behavior. The effective string tension 
<7 e ff is 10-20% reduced, compared to the string tension in the static case, in 
R < 0.8 fm and 0.5 GeV < M q < 1 GeV. The light-quark mass dependence 
of the reduction is also investigated. The effective string tension means the 
inter-two-quark confining force in baryons, and its reduction by the finite-mass 
valence quark effect is conjectured to be a general property for baryons and 
multi-quark hadrons. The significant change of the fundamental inter-quark 
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force is important not only for QCD but also for the quark-hadron and nuclear 
physics. 
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